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^ ABSTRACT 

We present a study of the possible plasmon excitations that can occur in systems where strong superconductivity is 
. present. In these systems the plasmon energy is comparable to, or smaller than the pairing gap. As a prototype of 
these systems we consider the proton component of neutron star matter just below the crust when electron screening 
is not taken into account. For the realistic case we consider in detail the different aspects of the elementary excitations 
when the proton, electron components are considered within the Random Phase Approximation generalized to the 
, superfluid case, while the influence of the neutron component is considered only at qualitative level. Electron screening 
plays a major role in modifying the proton spectrum and spectral function. At the same time the electron plasmon 
is strongly modified and damped by the indirect coupling with the superfluid proton component, even at moderately 
low values of the gap. The excitation spectrum shows the interplay of the different components and their relevance for 
each excitation mode. The results are relevant for neutrino physics and thermodynamical processes in neutron stars. 
O ■ If electron screening is neglected, the spectral properties of the proton component show some resemblance with the 
^ ' physical situation in high T^ superconductors, and we briefly discuss similarities and differences in this connection. 
^ , In a general prospect, the results of the study emphasize the role of Coulomb interaction in strong superconductors. 

PACS : 21.65.+f , 24.10. Cn , 26.60. +c , 03.75. Ss 
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lO ■ I- INTRODUCTION 

m 

The excitation spectrum in superconductors is strongly affected by the long range Coulomb interaction. In the sim- 
plest case of superconductivity in metals, the long range Coulomb interaction is responsible of the "Higgs mechanism" 
that produces a non-zero energy excitation at zero momentum. The otherwise present Goldstone mode, caused by 
the breaking of gauge invariance in a superconductor, actually disappears. Instead, a plasmon-like mode appears well 
above the gap energy, mainly at the plasmon energy of the normal metal [l', . In this way, the sound-like spectrum, 
below twice the energy gap, typical of e.g. neutral s-wave pairing superfluids, is modified to a spectrum that starts 
at finite high energy at zero momentum and varies smoothly with momentum. 

More complex is the situation in a high Tc superconductor, since electrons have no Fermi liquid behavior. However, 
. it has been argued 3] that along the so called c-axis plasmon-like excitations (" Josephsonplasmons" ) should be 
present, and indeed some authors have claimed that they can be experimentally observed [J, |5|. In this case, the 
plasmon frequency is smaller than the energy gap. 

The system over which we will focus our study is the homogeneous matter in neutron stars, just below the crust, that 
is expected to be superfiuid. In particular, homogeneous matter below Neutron Star (NS) crust is expected to have 
a proton superfluid component. The elementary excitations of the matter affect the whole thermodynamics and long 
term evolution of the star. Since the main components of the matter are neutrons, protons, electrons and muons @, the 
spectral properties of these excitations can have a complex structure. Collective modes in asymmetric nuclear matter 
have been studied previously, e.g. in Refs. 04^- In the astrophysical context, a study of the collective excitations in 
normal neutron star matter on the basis of the relativistic mean field method has been presented in Ref. llOl . The 
spectral functions of the different components in normal neutron star matter have been calculated in Ref. on the 
basis of non-rclativistic Random-Phase- Approximation (RPA) for the nucleonic components and relativistic RPA for 
the leptonic components. Different models for the nuclear effective interaction have been considered and a detailed 
comparison was done between some Skyrme forces and a microscopically derived interaction. 

In Ref. particular attention was payed to the plasmon excitations of the proton component, where the Coulomb 
interaction, among protons and between protons and electrons, plays a major role. In this work we will extend this 
study to the case of superfluid matter. The elementary excitations in superfluid neutron star matter have been studied 
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by several authors, often with controversial results |13l - ll8l |. It is well known that a neutral superfluid must present 
a Goldstone mode at low momentum, and this can have a strong influence on e.g. neutrino emission [Tst - fioj or mean 
free path. As already mentioned, for a single charged superfluid the Coulomb interaction suppresses the mode, which 
is replaced by the plasmon mode. In neutron star matter the physical situation is complicated by the multi-component 
structure, since the plasmon mode is in fact mainly an electron excitation, and furthermore the nuclear interaction 
couples neutron and proton excitation modes. We formulate the general theoretical scheme within the conserving 
approximations [2l|, which guarantee current conservation and the fulfilment of the related Generalized Ward 
Identities (GWI). However, we will focus the study only on the charged components, protons and electrons, while the 
complete treatment will be presented elsewhere in future works. The reason of this choice is that the role of Coulomb 
interaction is quite crucial in this case, while the coupling with neutron is expected to be weak [ll| . What is left out 
from the treatment is the possible role of the entrainment (2^ [23j when both neutron and proton components are 
simultaneously superfluid. However, below the crust neutron superfluidity disappears rapidly at increasing density, 
and we will therefore assume that only protons are superfluid. 



II. CONSERVING APPROXIMATIONS FOR SUPERCONDUCTORS. 

Let us remind first the conserving approximation scheme for normal systems as developed in Refs. (20l. [2]|. Since 
we are looking for collective excitations corresponding to density fluctuations, we are interested in the vertex function 
A, defined by 

A(12;3) =<T(p'(3)V/(l)^(2) > (1) 

where each symbol i = {vi,ai,Ti,ti) stands for coordinate, spin, isospin and time variables of a single particle state; 

are the creation and annihilation operators for the considered particles. For given spin and isospin values, 
the particle density is defined as p{i) = ■ip^{i)^p{i) and p'{i) = p{i) — < p{i) > is the density fluctuation. Following 
the standard conventions, the single particle Green's function G(12) is dcflned as 

G(12) ^ -t< r{7/i(l)^t(2)} > (2) 

It satisfies the Dyson's equation 

G-i(12) = 0^7^(12) -t/(12)-E(12) (3) 

where S(12) is the self-energy and Go the Green's function for non-interacting particles. We have also introduced 
a possible external single particle potential U, eventually local in space and time and acting on the particle density. 
The vertex function A of Eq. ([T]) describes the linear response of the system to such a local external potential 

where the functional derivative is taken at t/ = 0. By using the chain properties of the functional derivative, one can 
derive the integral equation for the vertex function 

A(12;3) = Ao(12;3) + Ao(12;P2')V(P2';45)A(45;3) (5) 

where a bar over a symbol i indicates integration and summation over the corresponding set of variables, and 

Ao(12;3) = -G(13)G(32) (6) 
The effective and irreducible interaction V can be expressed in terms of the functional derivative 

^(-^«' - -ii <') 

here the functional derivative is meant performed considering the bf Green's function as a functional of the external 
potential. 

If the self-energy S is approximated with a functional of the Green's function and of the interaction, both Green's 
function and self-energy must be calculated by self-consistent procedure. This is the first condition that a conserving 
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approximation must fulfil. In addition to the self-consistency, one has to impose a symmetry condition on the self- 
energy, that is automatically satisfied if it can be obtained as a functional derivative 

Here $ is formally a functional of the Green's function and of the interaction. Approximations of this type for the 
self-energy are called " ^-derivable" approximations and are automatically conserving [21]. It follows from Eqs. (|3l5l7p 
that a conserving approximation for the Green's function entails a conserving approximation for the vertex function 
A, provided the irreducible interaction V is obtained from Eq. ([7]). 

To make clear what a conserving approximation actually means, let us remind that the average density < n(r, t) > 
and current < j(r,i) > can be obtained from the single particle Green's function 

<n{r,t)> = Y.^ <ilj\r,a,t)^j{r,(T,t) >= iY.^G{x,x+) (9) 

<j(r,i)> = lE,A(Vr-Vr') <^^(r,(T,t)^/.(r',a,t) > |r'=r = E,(Vr- Vr')G(a;,a;+)|r'=r 
2m 

where = (r, a,t + e) (the isospin variable has been neglected), with e a positive infinitesimal quantity that enables 
to fix the correct time ordering of the operators in the equal time limit, according to the Green's function definition. 
A conserving approximation is such that the so obtained density and current satisfy the local conservation law 

^ <n(r,t) > +V <j(r,i) (10) 

The main reason of this result is the fact that, under the ^-derivable conditions, the contribution of the interaction 
term to the divergence of the current vanishes, and the conservation law of Eq. (jlOp then follows trivially from the 
kinetic term contribution. It has to be noticed that the interaction must be local, otherwise the conservation law 
cannot be written in such a local form. The Hartree and Hartree-Fock approximations for the self-energy are noticeable 
examples of conserving and ^-derivable approximations. These approximations generate the so called Random-Phase- 
Approximation (RPA) for the vertex function in the integral equation ([S]). Then this approximation for A conserves 
the current. 

The conserving approximation scheme was originally developed [2TI| for a normal system and for central particle- 
particle interaction. The theory can be readily generalized to a more complex interaction and to the case of superfiuid 
systems. 

First we have to extend the general formalism to the superfiuid case. The inclusion of superfluidity can be formally 
achieved by including in the above scheme a further discrete variable a labeling the destruction {a — 1) and creation 
(a = —1) operators, so that now the collective index is i = {ri,ai,Ti,ti,ai) = {xi,ai). Some care must be used 
in generalizing several relationships from the normal to the superfiuid case. In particular, the equation defining the 
inverse Green's function G(12)^^ must be written 

G-i(l;x3,a3)G(x3,-a3;2) - <5(1 - 2) (11) 

Then the single particle self-energy S(12) contains a normal S" and anomalous part E", corresponding to ai 7^ a2 
and ai = 02, respectively, and the Dyson's equation ([3]) still holds including the additional variable. Similarly, the 
generalized vertex function can be introduced by functional differentiation 

A(12;x3,a3,4) = ^TTtP^^^^^, = < r(^(l)V'(2)n '(^3, ag, a^,) > (12) 

where now n(a;3,a3,a3) = V (2^3, 0^3) (2^3, 013) is the generalized density, n' — n~ < n> the corresponding fiuctu- 
ation, and U couples linearly to n. Using again the chain property of the functional derivative one gets the integral 
equation for the generalized vertex function 



A(12;3) = Ao(12;3) + Ao(12;l'2')V(i:,2:,;45)A(45;3) (13) 

where a subscript ^ indicates that a sign change of the variable a has to be considered, e.g. 1* ~ (ri, cti, ti, — ai). The 
function Aq has the same expression as in Eq. ([6|) , provided the bf Green's functions are generalized to the superfiuid 
case. The correlation function A(12; 3) and the integral equation are the generalization to the superfiuid case 
of Eqs. ([1} and ([S]), respectively. Also within this extension of the formalism, an approximation is conserving if it is 
^-derivable and the Green's function is calculated in a self-consistent way. In the mean field approximation, which 
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generates the (generalized) RPA approximation for the vertex function, this means that the self-energy appearing in 
the RPA equation for the response function must be calculated self-consistently within the mean field approximation. 
For the pairing part this means that we have to use the BCS bf Green's functions. 

To illustrate this conserving approximation, we will check the conservation of current in the case of a pure pairing 
interaction . Let us take a zero-range interaction, that is the one we are going to use in this work 

?7pair = C/pair J dV^A^ (r , a) (r, -a) V(r , -a' ) V(r , a') (- 1) ^ (14) 

where spin-zero pairing has been assumed and Uq is the corresponding strength. The relevant component of the 
Green's function G{11') corresponds to ai = 1 and a'l = —1. This component will be indicated by G{xi,x[). The 
equation of motion for G can be written 

ih-^^Gixi,x[) + ^VlG{x^,x[) - S{t-t')5{r-r') (15) 

where xi — (r, (Ti,t) and x[ — (r',CTi,<'). The analogous equation for the derivative on the time t' reads 

^n-^G{x^,x[)-^WlG{xux[) = S{t~t')S{r-r') (16) 

+ zt/pairS. < T{^{r,ai,t)^\r' ,<j,tyj\r' ,-a,t')^{r' ,-<Ji,t')) > (_l)(i+-i-) 



Summing up the two equations of motion, performing the summation over the spin ai , and finally putting t — t' and 
r = r', the first two terms on the right hand side give the conservation law (Eq. (fTU)) '). because they correspond to the 
kinetic term contribution. The interaction contribution on the right hand side, that involves the two-body Green's 
function G2, factorizes in the mean field approximation. For the G2 contribution in Eq. (|15p one has: 

< T{i^\r,~ai,t)^{T,^cT\t)ij[v,a' ,t)^\r' ,ai,t') > « < T{'4j\r,-(Ji,t)^\v' ,ai,t')) >< ^{r,^a' ,t)ij{r,a' ,t) > 

- < (r, -CTi , t)V(r, a',t)>< r(V'(r, -a', t)ip^ (r', ai , t') > 

(17) 



A completely analogous factorization applies to the G2 contribution in Eq. (|16p. Taking into account the translational 
and time reversal symmetries 

<T{',p^{r,^ai,t)^j\r',ai,t')>= - < T{^{r,-ai,t)^{r' ,ai,t') > (18) 

the two factorized G2 are equal, but they are multiplied by an opposite phase in the equations of motion and therefore 
they cancel out. This conservation property is translated to the vertex function if the effective interaction V is 
calculated by the functional derivative of equation ([7]), properly generalized to include superconductivity. 

This property of the RPA equation, when only pairing interaction appears, was recognized a long time ago 01 ■ 
If other interactions are present, besides the pairing one, noticeably the Coulomb interaction, the argument can be 
repeated and this RPA property still holds, provided the additional interaction is local. The mean field correction to 
the self-energy in the mean field approximation amounts to a shift in the chemical potential, if the exchange term 
is neglected. The latter approximation is equivalent to the Hartree approximation, which is also conserving. If the 
interaction is non-local, or we include the exchange term, the local continuity equation (jlOp does not hold any-more. It 
can be restored in the effective mass approximation, where the effective mass is momentum and energy independent. 

It has to be noticed that the conservation law of the vertex function A, appearing in Eq. (I13p . is equivalent to the 
generalized Ward identities [2j . 

A direct consequence of the conservation law is the behavior of the density-current polarization tensor, defined as 

U,f,{x, x') = < T{j,ix)jf,{x')) > (19) 

where v, fi = 0, 1, 2, 3 and is the current four vector, with jo{x) = p{x) the particle density. This is the response 
function to a local external probe, which is linearly coupled to j^, provided the particle- hole interaction is local and 
the pairing interaction is zero range. It can be written in terms of the vertex function in a straightforward way 
The conservation law for 11 in the energy-momentum (w,q) representation can be written 



^ 11(0^) + ^1=1, 3 9i 11(1^) = 



(20) 
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For an isotropic system the second (current) term is proportional to for small |q|, and thus for any value of the 
energy lu the density-density component of the polarization tensor, in the limit of small |q|, vanishes at least as q^. 

To conclude this section, we will write explicitly the self-energy in the mean field approximation and the corre- 
sponding effective interaction for the RPA equation. Following the factorization of the two-body Green's function G2 
illustrated above, one can extract the self-energy for a general interaction v 



S(12) = 



/ [2Sa2,^ai < xia;i'|w|a;2a;2' >a {G{x2' ,tf , -1; xy ,l)5ai-i - G{xy ,tf , ~1; X2' ,l)5ai.i) 

+ <^a2.ai < XiX2\v\xvX2' >A G{xy , Ui, X2' ,t'^ , OLi)] (21) 



The effective interaction V, to be used in RPA equations, can be derived from Eq. ((2T|) : 
(5Xj(12) 

V(12;45) = ^^^^^^ i5{t2 - ti)5a2-ai5ai-i5ar,,i [< xiX5\v\x2Xi >A ^Qi.-i- < a::ia;4|w|a;2a;5 >a 5ai.i] 

1 

< X\X2\v\x\X^ >A (22) 

Again, if we assume that the interaction is the sum of a local particle-hole interaction and a local zero-range pair- 
ing interaction, the expression simplifies, and, as we will see, the energy-momentum representation is particularly 
convenient. 



III. EQUATIONS FOR THE RESPONSE FUNCTION 

In the BCS framework, the single-particle Green's function reads: 

G^^(fc) = a\\6a^-a2^„^^a^Gi}i)+a\6a^^a.2^„^-„^F{k)\ 



with 
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For simplicity we assume the pairing to be in the s-wave. We also consider the pairing interaction to be constant with 
a cutoff in momentum space, which corresponds to a contact interaction in coordinate representation. Furthermore, 
the interaction in particle- hole channel is local (it is a central density-density Coulomb interaction). In such conditions, 
the general equation ()13p for the vertex function reduces to an equation directly for the polarization function n(q, w), 
in its generalized form, appropriate to superconductors. As usual, the momentum representation is obtained by 
Fourier transform. From Eq. (jl3p one gets a system of four coupled algebraic equations to be solved for n(q, w); 
however it turns out that one equation is a linear combination of the others and the coupled components of the 
polarization tensor to be calculated are reduced to three. Finally, the generalized RPA equations read: 
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where the index S specifies that scalar (zero total spin) excitations are considered, and we have introduced the 
notation: 



X 



ph 



1 



X^p = XoFiq) ± XoFi-q) 
Xggc = l[X'J^i-q)-X^aGil)] 
The quantities X are the following four-dimensional integrals 

dk 



ph 
GG 



XPP 



(.q) 

gg(9) 



X'Jai-q) 

XGF{q) 
XGF{~q) 
XFFiq) 



(2^)4 
rlk 

dk 



G{k)G{k + q) ; 



G{k)G{-k-q) 



(2^)4 
dk 

^Mk)nk + q) 
dk 

J^,G{k)F{k~q) 
dk 

^^F{k)F{k + q) 



XpFi-q) = XpFiq) 



(24) 

(25) 
(26) 

(27) 



(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



They can be considered the different components of a generalized Lindhard function, suitable for the study of superfluid 
systems. Their explicit expressions are given in the Appendix. If the alpha indices are explicitly indicated, the 
polarization tensor can be written 



n(q,w) = n„i„'^.a2Q^(q,w) 

and the three components of the polarization tensor appearing in Eq. (j23p are defined as 



n(+) = -{lilUo.p+Ii-l-l-af}) 



(34) 

(35) 

(36) 
(37) 



where the values of the variables a and /3 are generic and are chosen according to the polarization tensor components 
that have to be calculated. 

In the so-called constant level density approximation, i.e. exact particle-hole symmetry, the off-diagonal matrix 
elements of Eq. (P^ . coupling the first equation with the other two, are of order q^. As a result, to order the 
equation involving the component 11^ is decoupled from the other two, and actually can be neglected. 

In the case of neutron star matter, one must take into account the coupling of protons and electrons by the Coulomb 
interaction. In this case the complete RPA equations form a system of three coupled equations for the polarization 
tensor, involving both proton and electron components 
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where X^ is the relativistic Lindhard function for electrons, that will be calculated in the Vlasov limit, and n'^^^^ is 
the corresponding part of the polarization tensor involving the electron (density) component. 
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IV. THE EXCITATION SPECTRUM 

The simplest situation is when only the pairing interaction is included. For the considered case of Eq. (j38p this is 
formally equivalent to put v^. — 0, and it corresponds e.g. to an uncharged superfluid. The corresponding expression 
for n^''^ can be easily obtained by solving the remaining two by two algebraic system. In the notation introduced in 
the previous section, and using the expressions of and Tl^^g^ given in Appendix, it reads 

/(I - C/pair) (39) 

One can verify that the numerator of this expression vanishes at q = for any non-zero value of uj, which implies that 
the response function , and the corresponding strength function, are proportional to for small q for all non-zero 
value of w. As shown in the previous section, this is a consequence of the conserving approximation we are following, 
which guarantees the conservation of current. Indeed, the continuity equation implies this property of the density- 
density component of the polarization tensor. This RPA property of current conservation is a well known result [4I 
in the theory of superconductors, where only pairing interaction is present, and it is also a consequence of gauge 
invariance. 

According to a general theorem on broken symmetry, the strength function for uj < 2A is characterized by the 
presence of a "Goldstone boson" , i.e. a phonon-like excitation with an energy proportional to |q| for small enough 
|q|. It is well known Q that in this limit the velocity of this mode is vp/y/3, provided the weak coupling approximation 
is valid, where vp is the Fermi velocity. The energy of the Goldstone mode is obtained just by the equation 

1 - C/pair = (40) 

which corresponds to the vanishing of the determinant of the matrix on the left hand side of Eq. ([38)) . As expected, 
according to Eq. ([39]), the strength function has a delta singularity at the energy of the mode, since it is undamped 
in the considered limit of pairing interaction only. 

It is a classical result [1] that in a charged superconductor the Goldstone mode cannot be present since the standard 
theorem on broken symmetry (the gauge invariance in the present case) is not any more valid if the particles interact 
also by a long range interaction (the Coulomb potential). The mode is replaced by the plasmon mode. The plasmon 
excitation has an energy which starts at a non zero value at q = (the plasmon energy) and varies slowly with the 
momentum. This conclusion can be obtained by considering the first two equations psp . including Vc and neglecting 
the third equation and the electron component. We deal then with a charged superconductor. One finds that the 
determinant of the matrix on the left hand side of Eq. p8|) in this limit vanishes at non-zero energy for q = 0, since 
the factor in front of actually vanishes, as it happens in a normal charged fluid. 

Here we focus on the case when the pairing gap is comparable to or larger than the plasmon frequency, as it can 
happen in high Tc superconductor. It can be true also in neutron star matter for the proton component. Even 
if in the latter case the electron screening changes the physical situation, we will discuss the neutron star proton 
superconductivity, stressing the possible analogy with the high Tc superconductors, where the physical situation is 
more complex due to the crystal electronic structure. Then we will introduce the electron screening and discuss the 
effect of proton superfluidity on the electron plasmon. 

The first main point that we want to address in a general framework is where the position of the plasmon excitation 
can be expected at increasing value of the pairing gap. At first sight one could expect that the original excitation 
spectrum, present in the normal system, is just shifted upwards by an amount comparable with the energy gap in 
the superfluid system. However it turns out that this is not the case. Let us consider the neutron star matter at a 
total baryon density equal to the saturation one, i.e. 0.16 fm~'^. At this density the proton fraction is expected to be 
3.7% [Ul, l12], and the proton plasmon energy Up « 2.105 MeV. Under these physical conditions, the full structure of 
the excitation spectrum is illustrated in Fig. [l] for different values of the pairing gap. 

The branches of the excitation spectrum are defined as the zeroes of the determinant of the real part of the 
matrix at the left hand side of Eq. psp in the energy-momentum space (w, |q|). The density and gap values can be 
considered realistic [Tlj . [l2l . [23 |. Some comments are in order. When the gap A is much smaller than the plasmon 
energy, the structure of the excitation branches resembles the " thumb like" shape typical of the charged normal Fermi 
liquid [lll,[ll [11-113]. In this case the lower branch is over-damped, and indeed it does not correspond to an actual 
excitation, while the upper one is just the plasmon mode, that is mainly undamped up to a momentum nearby the end 
of the "thumb", above which no excitation mode exists. In the considered superfluid case, however, the lower branch 
must start above the forbidden energy region ui < 2A. At increasing value of A, the spectrum changes considerably. 
Not only the momentum dependence is modified and quite different from a quadratic shape, but the whole spectrum is 
compressed and stays just above 2A. When 2A approaches the plasmon energy ujp of the normal system the spectrum 
can exist only in a narrow range of momentum and finally no excitation is possible. Therefore no plasmon excitation 
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FIG. 1: Proton excitation spectrum at different values of the pairing gap, as indicated by the labels. Electron screening is not 
included. 



is possible inside the forbidden region, at variance with high Tc superconductors, where a plasmon mode is observed 
well below 2A along the c-axis. However in high superconductors the reason of the low value of the plasmon mode 
is due not only to the low density of the charge carriers but mainly to the very high value of the insulating dielectric 
constant along the c-axis [1, Q]. The coherent superfluid motion is possible because the pairing coherence length 
is longer than the inter-plane distance in the c-axis direction, so that the planes form a micro-array of Josephson's 
junctions (hence the name of " Josephson plasmon"). 

Going back to the neutron star matter, we notice that if the electron component is introduced, the overall structure 
of the spectrum is modified. In fact, electrons are much faster than protons and are able to screen the proton-proton 
Coulomb interaction at all considered frequencies. Since then the proton-proton effective interaction is of finite range, 
a sound-like branch appears again below 2A. This excitation can be considered a pseudo-Goldstone mode. It is 
determined by both pairing and screened Coulomb interaction. In fact, without the pairing interaction a sound mode 
is anyhow present at low momenta IjJ; it is strongly damped, and its velocity is determined by the screened Coulomb 
interaction. As discussed above, if only the pairing interaction is considered, a Goldstone mode is necessarily present 
below 2A, essentially undamped, with a velocity that in the weak coupling limit is independent from the interaction 
strength and equal to wf/\/3. The pseudo-Goldstone mode turns out to have a much higher velocity, about three 
times larger. 

The structure of the excitation spectrum is illustrated in Fig. [2] for A = 1 MeV and at the same total baryon 
density. The reported branches correspond to the ones where the proton strength is dominant. One notices the 
pseudo-Goldstone mode below 2A and the pair-breaking mode close to 2A. At increasing momentum the pseudo- 
Goldstone branch deviates from the linear behavior and approaches slowly 2A and the two branches seem finally to 
touch at an end point of the spectrum. 

For much smaller value of the gap the spectrum looks slightly different. The pseudo-Goldstone mode and the 
pair-breaking mode actually cross. Above the crossing the upper branch corresponds to the sound mode present of 
the normal system, while the lower branch corresponds to the over-damped mode that is also present in the normal 
system. Roughly speaking, the "thumb-like" shape of the spectrum above 2 A resembles the similar shape that is 
present in the normal system [ll], [T^ . We will not discuss in detail the evolution of the spectrum as a function of the 
gap value since this is outside the scope of the present work. In any case, if the value of the gap is vanishingly small, 
one recovers the spectral shape of the normal system. 

Up to now we have neglected the nuclear interaction. The inclusion of the proton-proton and proton-neutron 
effective interactions in the description of the collective modes in superfluid neutron star matter is beyond the scope 
of the present work, and we will keep the discussion at a qualitative level. In previous papers [ill [l^ we have shown 
that the neutron-proton coupling in the very asymmetric matter of neutron stars is weak, and therefore the proton 
excitations are approximately decoupled from the neutron component. The nuclear proton-proton interaction can of 
course modify the velocity of the pseudo-Goldstone mode, but the above scheme of the excitation spectrum should 
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FIG. 2: The branches of the proton excitation spectrum for A 



1 MeV. Electron screening is included. 



be still valid. The detailed study of the general case of superfluid neutron starmatter is left to future works. 



The spectrum of the electron component in neutron star matter is characterized by a plasmon excitation. Since the 
electron Fermi velocity (very close to the speed of hght) is much higher than the proton Fermi velocity, the protons 
act essentially as a static positive background. However some dynamical electron-proton Coulomb coupling is still 
present and produces a relatively small shift of the electron plasmon energy with respect to the purely static limit. At 
the plasmon energy some proton strength is then present, but it turns out to be in general quite small [ill. [l2j|. If the 
protons are superfluid another phenomenon occurs, that is not present in the normal system. The electron plasmon 
is coupled by the Coulomb interaction to the particle-hole proton excitations, characterized mainly by a sound mode. 
If the sound mode energy is higher than 2A, it can decay since it is in turn coupled to the pair-breaking mode by 
the pairing interaction. This indirect coupling of the electron plasmon to the proton pair-breaking mode produces 
a damping of the excitation, and therefore its width increases at increasing value of the pairing gap. This can be 
studied by considering the electron strength function 



It gives the strength of the electron component at a given excitation energy E and momentum |q|. 

In Fig. [3] is reported the electron strength function for |q| =5 MeV as a function of E and for different values of the 
pairing gap, again at the total baryon density p = 0.16 fm~^ [Til [T^. One can notice the sharp plasmon peak that 
is obtained for the lowest value of the pairing gap, as expected for a normal system. However, even for moderately 
small values of the gap the width of the plasmon increases rapidly and its strength is spread in a larger and larger 
energy interval, corresponding to a stronger damping. In fact the energy-weighted sum rule for the electron strength 
must be satisfied, as we checked numerically, and therefore the electron strength is also present at much lower and 
higher energy with respect to the plasmon energy in the normal system. This can have some consequence for e.g. the 
neutrino emission in the cooling process of neutron stars. 



V. 



THE ELECTRON PLASMON IN NEUTRON STAR MATTER. 




VI. SUMMARY AND OUTLOOK 



We have studied the possible plasmon excitations in strong superconductors, where the pairing gap A is comparable 
with the plasmon energy. As a study case we have considered the proton component in the homogeneous neutron star 



10 



0.008 



0.006 



0.03 


1 ■ 1 ■ 


n riR 




0.12 








J 


l>- 



6 7 8 9 

E (MeV) 

FIG. 3: Electron strength function at different values of the pairing gap, as indicated by the labels. 



matter. We found that the plasmon excitation disappears if 2A is equal or larger than the value of the plasmon energy 
in the normal system. This should be a general result, not connected with the particular features of the neutron star 
matter. This finding is in contrast with what is predicted and observed in high-Tc superconductor, where a plasmon 
excitation is predicted and observed well below 2A. However, the physical situation is quite different in these solid 
state systems, where the smallness of the plasmon energy is not due to the low density of the charge carriers, but to 
the high value of the insulating dielectric constant. 

In a more realistic description of the neutron star matter one has to include the electron component. This modifies 
the spectrum, since the electrons are faster than protons and screen the Coulomb proton-proton interaction. The 
plasmon is replaced by a pseudo-Goldstone excitation below 2A. Finally, the inclusion of the particle-hole nuclear ef- 
fective interaction among protons and neutrons is expected to further modify the spectrum, but some of the qualitative 
features of the proton component excitations are expected to remain valid. 

We have then considered the electron plasmon excitation, that in neutron star matter occurs at energy much higher 
than 2A. Although the proton-proton pairing interaction can affect only indirectly the electron excitations, it turns 
out that even at moderately small value of the pairing gap the electron plasmon is strongly damped. 

On general ground, our results stress the relevance of the Coulomb interaction in strong superconductors, in 
particular in neutron star matter. 



Appendix A: Generalized Lindhard function and numerical method 



The relevant functions X appearing in the kernel of the RPA equations, reported in Eqs. (|24II26|) . after the integra- 
tions over energy, can be written 
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Despite the approximation of a constant pairing gap, these integrals cannot be done analytically and their numerical 
evaluation requires particular care. We follow the method of assigning to the quantity rj a small value, simulating 
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an infinitesimal, and then integrating by an adaptive method the remaining two-dimensional integrals. We used the 
subroutine DT20DQ of the Visual Fortran package. In this way both real and imaginary part can be calculated in 
few minutes in a simple PC. We found that the results are quite stable if we take for r] values between 10~^ and 10~^ 
(energies are all calculated in MeV). For production calculations we used the value 1.5 x 10~^. 

Finally, for completeness we give here the expression of the free polarization functions appearing on the right hand 
side of RPA Eq. ((38|). in the case we are interested in the density-density response function, i.e. the spectral function 
( a -1 and ^ 1 ) 
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Using this expression, for the pure pairing case (i.e. Vc — ^ 0) one gets the expression of Eq. p9l) for the density-density 
component 11^''''. Then, putting q = in the previous expressions for the functions X, one finds 
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(27r)3 E^[{2Ek) 



(A7) 



The relation shows that Ilg''^\q,qo) vanishes at q = for any value of qo, i.e. it is proportional to q^ in the small 
|q| limit for all qq. This result can be readily generalized to the case where a finite range particle- hole interaction is 
present. 
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